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Abstract. We present a minimal statistical model designed for the description of rare- 
hadron multiplicities in nucleus-nucleus collisions at energies below the threshold of the particle 
production in binary elementary collisions. Differences to more conventional canonical statistical 
model are explained. The minimal statistical model is applied to the description of multiplicity 
ratios involving H hyperons, which are measured by the HADES collaboration at GSI-SIS. It 
is argued that the HADES data cannot be reproduced by the model based on the statistical 
equilibrium and the strangeness conservation. The data remain underpredicted even when in¬ 
medium potentials acting on hadrons are taken into account. This hints to non-equilibrium 
production of the H hyperons and their continuous freeze-out. 


1. Introduction 

Data obtained by HADES collaboration pQ show that doubly strange H hyperon can be 
produced in nuclear collisions at energies well below the production threshold for nucleon- 
nucleon collisions. In Ar-|-KCl collisions at Ebeam = 1-76 AGeV the observed ratio of H to 
(A -|- T,^) multiplicities is AIh/AIa-i-so = (5-6 ± x 10“^, and some H’s are also seen in 

Au-|-Au collisions at 1.23 AGeV. Statistical model with canonical suppression of strangeness 
production predicts a value of the order (10“^) for the ratio. 

2. Isospin symmetry and total strangeness content 

The colliding nuclei with the atomic number A have usually different number of protons (Z) and 
neutrons (A — Z) and thereby the matter created in their collisions has a non-zero total isospin 
projections quantified by the asymmetry coefficient rj = {A — Z)IZ being different from one. 
In our calculations we assume that this asymmetry is reflected in the observed multiplicities of 
hadrons belonging to one isospin multiplet, e.g, 
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For the Ar-|-KCl collisions we have r/ ~ 1.14. 

We can use the relations 0 to determine the total amount of produced ss quark pairs 
(distributed among hadrons). In a baryon-rich environment, created in collisions at SIS energies, 
the s antiquarks can only be accommodated in kaons. Furthermore, under such conditions kaons 
do not find partners for inelastic scattering and once produced they are guaranteed to leave the 
system. Thus the number of ss pairs can be reconstructed as 

Mss = Mk+ + Mko = (1 + 'n)-^K+ ■ 


( 2 ) 






Once the numbers of produced s and s quarks are known we assume in our minimal statistical 
model [2] that the s quarks are distributed among strange hadrons (anti-kaons and hyperons) 
with relative abundances determined by thermal and chemical equilibria and by the strangeness 
conservation. For our analysis we use the ratios which formally do not depend on the produced 
strangeness [2] 
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These ratios has the same total number of strange an anti-strange quarks in both numerator 
and denominator. 


3. The minimal statistical model 

At subthreshold energies no strange particles are created in incident nucleon-nucleon collisions 
and all strangeness must be produced during the evolution of the hot hreball created in the 
collision. Creation of strangeness is a very rare event, with multiplicity 

Mk+ = (2.8 ± 0.4) X 10-2 . (4) 

The amount of ss pairs is thus far below the saturation and can be treated perturbatively. It 
will not only scale with the volume, but also we can assume that it will steadily grow with time. 
If — for simplicity — we assume ideal hydrodynamics, then the only length scale is given by 
where V is the fireball volume, and it must determine also the lifetime [3]. At fixed impact 
parameter, the probability of creating one ss pair is given as 

W = . (5) 


For a fixed volume V, the probability to have n of the ss pairs is Poissonian 
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The HADES data have been obtained after averaging over collisions with different centralities. 
The dependence on 1/^/3 ^ makes the averaging non-trivial. If we denote the centrality 

averaging by angle brackets, (•••), we can relate the parameter A to the measured kaon 
multiplicity as follows 

(IF) _ Mss _MkA^ + v) 
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In the minimal statistical model we assume that in an event with volume V and n pairs of 
ss quarks, the s quarks are distributed into hadrons according to the probability 
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This is the probability to hnd hadron a with the mass nia, baryon number Ba and strangeness 
content Sa- The chemical potential is determined by the fireball freeze-out density taken 0.6po, 
where po = 0.16 fm-^ is the nuclear saturation density. The probability Q is normalised — for 
each n separately — by the factor . 
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Figure 1. Multiplicity ratios as functions of the temperature. Red dotted lines with purple 
bands show the measured data values ([^. Dashed lines are minimal statistical model predictions 
according to (10). Solid lines are predictions with in-medium potentials included. For the 


M-k- jM.K+ ratio we show results with two values of K optical potential acting on U : —75 and 
— 150 MeV (other curves are calculated with —75 MeV). Dash-doubly-dotted curves show ratios 
with only nucleon potentials included. 


In total, average multiplicity of hadron species a is then given by the sum 

Ma = + (Mf)) + (Mf)) + ... (9) 

where the r.h.s. lists all contributions to the average multiplicity from events with one, two, 
- in) 

three, etc. ss pairs. Note that Ma is understood as the number of hadrons a produced in 
events with n ss pairs, divided by the number of all events. These multiplicities are determined 
with the help of Pa ' and Pgs- The H hyperons, caring two s quarks, are produced only in 
events with two and more ss pairs. 

4. The results 

Our result for the H multiplicity ratio is 

Ms _ {V^/^){V) Ps/{Pk+Pa+Pe) . 

In the pre-factor we collected the modifications with respect to the “usual” statistical model 
with the canonical suppression. This and all other calculated ratios are displayed in Fig. 
by dashed lines as functions of temperature. We see that the calculated Ms/Mj^^-soMx+ is 
clearly below the measured error band. In Fig. we also revised the data on ratio based 

on the actual isospin ratio from ([^, see [2]. 

By trying to improve the agreement of our calculations with the data we have introduced 
density-dependent scalar and vector potentials Sa and 14 [Sj- They modify the statistical 
distribution function by replacing the dispersion relation to Ea{p) = \/{nia + -|- p'^+Va- The 

chosen values were Vat = Fa = |Fa = |Fs = 3Fh = ISOMeVps/pO) andS'Ar = — IDOMeVps/po- 

The scalar potentials for strange particles follow from the relation Sa = {Ua — Va{po))pB/Po, 
where we use Ua = —27MeV, C/s = 24MeV, Us = —14MeV, C/^ = —75 or —150 MeV and 
= 0. This leads to changes in our theoretical results which allow to explain K~/K^, A/K~^ 
and Yj/K^ ratios by one freeze-out temperature ~ 70 MeV. However, the changes are not large 
enough to bring the H/A/it'"*' ratio into agreement with the measured data. 

It is interesting that accounting for the trigger on central collisions actually makes the 
disagreement even bigger. The culprit is in the volume factors which reduce the overall result 
when selection on larger volumes is made. 




























5. A prediction for Au+Au collisions at 1.23 GeV 

Although the model disagrees with the Ar+KCl data, for curiosity we can make a prediction also 
for the Au+Au collisions recently recorded by HADES. As seen from Fig. the temperature 
dependence of Ads/Af a_|_soA1x+ is very weak. The crucial difference to the smaller system is 
in the volume which is almost five times bigger for Au+Au than for Ar+KCl. As the ratio (10) 
is proportional to V~^, we expect that it will be five-times smaller in gold-induced reactions 
than in Ar+KCl. Below the strangeness production threshold this statement is independent of 
the collision energy. 


6. Conclusions 

In comparison to the canonical statistical model, the minimal statistical model assumes that 
strangeness must be produced during the hreball evolution and takes into account the proper 
scaling of rare species multiplicity with the volume. This leads to non-trivial differences 
between the two models and the result of the minimal statistical model on centrality-averaged 
multiplicities depends on a set of averaged volume factors. 

In any case, the conclusion we could draw from our study is that cascade production is non- 
statistical. We speculate that it might be caused by immediate decoupling of the produced 
S’s from the bulk matter. At an early hot stage some S’s are produced which are later not 
annihilated at lower temperature. Those particles chemically decouple from the system. 

We speculate that the most viable for such a production of H would be the strangeness 
recombination reactions, either involving K and a hyperon, or two hyperons. Note that we have 
provided parametrised cross-sections for such reactions in [1] . An interesting options appears in 
the KA channel, which in vacuum has a high reaction threshold of 1609 MeV. Due to effective 
lowering of kaon mass in the baryon-rich environment it can decrease so that eventually it hits 
the resonance around 1530 MeV. This could increase the cross section for S production 
considerably. 

Kinetic calculations aiming the the validation of these mechanisms are planned for the near 
future. 
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